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ABSTRACT 

In this paper, present computational methods for solving numerical multiple integrals with constants and variables 
limits of integrations which are based on Haar wavelet, Block-pulse function, Chebyshev wavelet and Sinc-function 
methods. This approach is the generalization and improvement of these methods and compared with Gauss Legendre 
quadrature methods. The main advantages of the generalized method are its more efficient and simple applicability than the 
previous methods. An absolute and relative errors are estimated of multiple integrations are considered up to three 
dimensions. Finally, we also give some numerical examples to compare with existing methods and the benefits of proposed 
methods have to find their computation efficiency. 
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1. INTRODUCTION 

Integrals and derivatives are fundamental to calculus. The integrals involve in mathematics, physics, and 
engineering applications. To find the solution of integrals is preferred clear procedures. The analytical solutions are not 
helpful for the real-time applications. Some integrals cannot be found with exact solutions analytically, if some special 
functions involving the integrals which is a challenge to compute, and difficult to finding the accurate solution is too 
time-consuming. This is the main reason to go for numerical methods for approximating integrals, nowadays to compute the 
integrations using digital electronic computers. So numerical integration is the best ways to find most approximate solutions 
of the integrations. In last decades, wavelets are developing tool in this area. There are different types of wavelets and 
approximating functions have been used in numerical approximations, some examples are the Haar wavelet [Maleknejad and 
Mirzaee (2005)][Aziz and Haq (2010)], Block-pulse function [Rabbani and Nouri (2010)], Chebyshev function [Yuanlu 
(2010)] [Rostami et al. (2012)] and Sine function [Stenger (2012)]. The solution of numerical integral has been done up to 
triple integral using Haar wavelet [Ahmedov and bin Abd Sathar (2013)] [Aziz and Haq (2010)]. The single integral has been 
done using Block-pulse function, Chebyshev function and Sine function [Rostami et al. (2012)]. 
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In this paper, we generalized as N dimension integral (N = 1,2,3 ,....) using the proposed methods, some examples 
of numerical multiple integrations (consider only single, double and triple integrals) with variable and constant limits are 
provided to show the precision of proposed methods and comparison between them. Those methods are compared with 
Gauss Legendre quadrature method. Absolute and relative errors are estimated of numerical multiple integrations. 
The advantage of proposed the methods has to find their computation efficiency. 

In associate with this paper as follows. In section 2 as Numerical technique for single integrals of proposed methods, 
and the subsections 2.1 to 2.4 as method of numerical integrations based on Haar wavelet [Ahmedov and bin Abd Sathar 
(2013)] [Aziz and Haq (2010)], Block-pulse function [Rostami et al. (2012)], Chebyshev wavelet, Sine function 
[Rostami et al. (2012)] respectively. In section 3, the method for generalized N dimension integrals. In section 4, has given 
numerical examples with results and section 5, is given results and discussion. Concluding remarks are given in the last 
section 6. 


2. NUMERICAL INTEGRATION FOR SINGLE INTEGRALS 

2.1. Haar Wavelet 


The Haar wavelet is a sequence of rescaled square-shaped functions which together form a wavelet family or basis, 
and its scaling function 0 x (x) can be described as on the interval [a,b) 



1, for x E [a, b ) 
0, otherwise 


( 1 ) 


The Haar wavelets mother wavelet function 0 2 0*0 can be described as on the interval [a, b) is given by 
f 1 for x E [a, 

4> 2 W = | _! forx e (2) 

V 0 otherwise 

In general, Haar wavelet family defined for x E [a, b ) except the scaling function can be described as 
( 1 for x E [a, (3 ) 

4> iC*0 = j-l for x E [/3, y) (3) 

v 0 otherwise 

for all i = 3,4,5,.,2 M where a = a + = a + ^—^(k + 0.5), y = a + ^—^(k + 1) and M = 2 J , J be 

m m m 

the maximum level of Haar resolutions and k = 1,2,.,2 M 


2.1.1. Method of Numerical Integrations Based on Haar Wavelet 

Any function /(x) which is square integrable function on the interval [a,b) can be expressed as an infinite sum of 
Haar wavelets 


fix') = Zi=i a i ( Pii x ) (4) 

In practice, only the first 2M terms of the sum is considered as, 

fix) = Yd=i a-iCpiix) (5) 

Lemma 2.1 The approximate value of the integral is 
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$a /( x )^ x = a i(b - a) (6) 

Proof. 

= o 

for all i = 2,3, . 

and 

fa 0i (x)dx = b — a 

£f(x)dx = Y.f=i a-i I* 4>i(X)dx 
= %(£> — a) 

In connection with the above equation Haar approximate only one co-efficient in the evaluation of definite integral 
to calculate the Haar co-efficient a l9 we consider as a nodal points as, 

= a + ^(fc + °.5) (7) 

the discretized form of eqn. (5) can be written as 

/(*k) = Yd=l a i<Pi( x k ) (8) 

Lemma 2.2 The solution of the system for a k is 

% = f(x k ) (9) 

Proof. We can prove by an induction method on J, where M = 2/ for / = 0, M = 1 

/Oh) = % + a 2 

f(x 2 ) = a 1 -a 2 
therefore 

ai=^[/0i) + /0 2 )] 

this is true for / = 0, we assume that the lemma is true for upto / = n — 1, let us consider J = n the new system 
has contain 2 n equation involving 2 n variable, and replacing 2 a k by a k and f(x 2k _ 1 ) + f(x 2k ) by g{x 2k ), thus we 
have 


a'k = 9(x k ) 

substituting a k and g{x 2k ) values 

a k = /(**) 

a k = ^Ifc=i f(x k ) 

the lemma is true for J = n, in particular 


( 10 ) 
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«1 =^Z&/(**) 


-Xk-i/(a + — (fc + 0.5)) 

2M y v 2M v " 


substituting equations (11) and (7) in eqn.(6), the formula for Haar wavelet is 

Sa f(x) d X 

2.2. Block-Pulse Function 

Definition 2.1 An m-set ofBPF is defined as follows 
fl,forte [<izlE, E) 

— J m m 


4>i(x) = . 

0, otherwise 
with t E [0 ,T), i = 1,2,3,. ,ra and h = — 

L m 

2.2.1. Method of Numerical Integrations Based on BPF 

rb 

We consider the integral J a f (pc)dx by using x = (b — a)t + a, we have 

J* f(x)dx = (b — a) g f((b - a)t + a)dt 
Theorem 2.3 The approximate value of the integral is 

gfwt~±zr =1 fi 

Proof the orthogonal property of BPF is the basis of expanding functions into their BPF series as 
fit) = Mm 
Jo f(t)dt = YT=i /Jo' 0t(Odt 

~ m Li=1 Ti 

Let us consider the nodal points as 

2k—1 


t k = ——,/e = 1,2,. ,m 

K 2m ’ 


fih) = 5T=1 Miih) = fk 

The approximate value 

/„■ fwt * iir-, 

In general 


L,f(t)dt 


^-ZT-1 fi a + — (i - 0.5)) 

m m 


if m = 2M, M = 2f the BPF becomes Haar wavelet (see Tables) 

2.3. Chebyshev Wavelet 

Definition 2.2 The Chebyshev wavelets , cp nm (x),n = 1,2, ., 2 fc_1 and m = 0,1,_ ,M — 1 


(ID 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


(17) 

is defined on the 
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interval [0,1) as, 


071,771 00 


where 


Tm(t) = 


2I? m (2 fe x -2n + 1), for x € 
( 0, otherwise 


—, for ra = 0 

V7T 


(18) 


JM* 


for m > 0 


and k is any any positive integer and m is the degree of Chebyshev polynomials of the first kind. Here, T m (t ) is the 
Chebyshev polynomials of the first kind of degree m which is orthogonal with respect to the weight function W (t) = == 
on the interval [-1,1] and satisfy the following recursive formula 

f 

2.3.1. Method of Numerical Integrations Based on Chebyshev Wavelet 


(T m + 1(0 = 2tr m (t) — T m _ 1 (t'),m = 1,2,3,. 
T 0 = landT^t) = t. 


An integral over [a, b] must be changed into an integral over [0, 1]. This change of interval can be done in the 
following way 

C f( x ) dx = (b-a) g f((b — a)t + a)dt (19) 

Theorem 2.4 The approximate value of the integral is 

i _k 

/o 1 mdx * nil [c„, 0 + n=i c n , 2l ] (20) 

Proof. Any function /(x), which is square integrable in the interval x E [0,1), can be expressed as 

/(X) — YiTl= 1 2777=0 ^77,777077,771 (x)X E \CL, b) 

so we can approximate /(x) as 

fix) = nil Zm =0 C n ,m<Pn,m( x ) x £ [a,b) 

fo f( x )dx = ZUl £m=0 C n,m f) <Pn,m( x ) dx (21) 

let us consider 

fo (pn,m( x )dx = 21 fH T m (2 k x -2n+ 1 )dx 

2 k-i 

2~f_\ T m (t)dt 


since 


i r o.n 

/_! T m (t)dt = ] -2 

(•nt 2 -1 ’ 


0 ,m is odd 
m is even 
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we have 


Iq 




^1 

O 

II 

s 

V7T 

0 

m is odd 

3 -k 

2 2 

m is even 




substituting the above equations in eqn. (21), then we get eqn. (20). 

Now to calculate the coefficient 0 ) and Q n 2 i) °f Chebyshev wavelets, and consider the nodal points as 


x = q = 12,., 2 K ~ 1 M 

q 2 k M ^ 


and 

f(Xq) = Hn=l £m =0 C n ,m4>n,m( X q ) 

From the definition of Chebyshev wavelet 


(n- 1) < 2q—l n 
2 k-i ~ 2 k M ^ 2 k ~ 1 

SO 


q = (n — 1)M + i i = 1,2,3, . ,M 


/(*«) = / 


and 


k 

0n,o(X?) = ~7E 


l+fc 


V5 M 

The above system of equations changed to the below system of equations is given as 


k 

22 


f . 2(n-l)M+2t-l . _ 22 ■ yM-l 2 T 

2 k M j +2,771=1 ^ l mK M 


l+k 

? m-i 2 2 r r 21 ' 1 
ViF ^ m 


for all i = 1,2,3, . ,M 

The coefficient C^ n 0 ) and C^ n 2 i) can be easily calculated from the above system of equation. 

2.4. Sine Function 


Definition 2.3 The Sine function is defined on the whole real line by 


sinc(x) = 


sinCroO for ^ ^ q 

TCX 

. 1, for x = 0 


For any h > 0, the translated Sine function with evenly spaced nodes are 
s(/, ft)(x) = sine J = 0, ±1, ±2,. 


( 22 ) 
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This is called ] th Sine functions. If f is defined on the real line, then h > 0. The series 

whenever the series f is converging. The expansion of f is called Whittaker cardinal expansion. 

2.4.1. Method of Numerical Integrations Based on Sine Function 

In this case we take the limits (a,b), where —oo<a<b<oo 


0(x) = Ing) 


0'O) = 


b-a 


the basis function on (a,b) is given by 
S(J, /i)</>(x) = sine Jh ) 
we consider the nodal point as Xy E (a, b) 


Xj = (p 1 (Jh) = 


a+be ]h 

l+eJ h 


( 23 ) 


(24) 

(25) 


(26) 


Theorem 2.5 Let L a (D) be the set of all analytic function, D be a simply connected domain in the complex plane, 
E L a (D) with 0 < a < 1 and 0 < d < n, let N be the positive integer, and h = (^)2 and there exist a positive 
constant C 2 which is independent ofN, such that 


C fV)dt - /l2; 


f(xj} 


< c 2 e (_Jld “ w)2 


I = ~ N *\xj) 

so the approximate integral based on Sine function is 




(27) 


3. METHOD FOR GENERALIZED N DIMENSION INTEGRALS 


Consider the multiple integral with variable and constant limits of the type 

/a nr = 2 f( Si , Sl )d Sl d Sl (28) 

with s t = s 2 ,s 3 , . ,s n 


Applying the above proposed methods of weights (w 7 ) and the corresponding nodes (x ; ) in the above equation to 
the inner integrals and by treating the outer integral variable s ± as constant, the above equation can be written as 


ffr pm r d ( s i) 
Hi=2 Jc(sj) 


f ( Si, s 1 )ds i ds 1 



Lj=i 


w j( s i)f ( X J> s i-i)ds 1 


let us consider 


a(si) = nr=2 zt= i w j( s i)f( x j’ s i- 1) 

again applying the above proposed methods of weights (w ; ) and the corresponding nodes (x ; ) in the above 
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equation to the outer integral leads to 


rrn r d ( s i) 
1a lh=2 Jc(Sj) 


f(s i ,s 1 )ds i ds 1 


fa 9 (s 1 )ds 1 ~ Yfh Wj,g(Xj) 


(29) 


The domain of integration for the above methods is conventionally taken as [0, 1) except Gauss Legendre 
quadrature method and this method is taken as [-1, 1]. Sum of weights is a range of domain, n be the number of dimensions 
and m be the number of points. In this approach, the same number of points is considering to each integrals. 


4. NUMERICAL EXAMPLES 


The following examples are given to show the accuracy and efficiency of numerical integration by using Haar 
wavelet, Block-pulse function, Chebyshev wavelet and Since function are compared with Gauss Legendre quadrature 
method and exact solution. The absolute and relative error has been calculated and compared to each method for the given 
examples are shown in the below tables and to illustrate the computation performance of multiple integrals in Figure 1 to 6 
obtained with above method. 

Definition 4.1 Suppose that x t and x a denote the true and approximate values of a integrals then the error 
incurred on approximating x t by x a is given by 


e = x t -x a 

i.e. magnitude of the absolute error e a is given by 

e a = \x t -x a \ (30) 


Definition 4.2 Relative Error or normalized error e r in representing a true value x t by an approximate value x a 
is defined by 


_ \x t -x a \ 
r \Xt\ 


(31) 


4.1. Examples for Single Integration 
Example 1 

fj*" Vxcos 2 (xz)dx 


The solution of the integral is 4.531824108052927 Using eqn. (30). we calculated the absolute error, where x a is 
the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and exact 
solution are shown in Table 1. and using eqn. (31). we calculated the relative error and shown in Table 2. From the table as 
considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 1. obtained with 
above method. 
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Figure 1: Time Evolution of Given Methods for Example 1 


Table 1: The Absolute Error for Example 1 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 


(m,NN) 


(k,M,NN) 


(N,NN) 

e a 

(n,NN) 

e a 

2,8 

4.0009E-01 

10,10 

1.8978E-01 

3,2,8 

4.0009E-01 

4,9 

8.4416E-01 

4,4 

6.4792E-03 

4,32 

1.4264E-02 

30,30 

1.6140E-02 

4,4,32 

1.0206E-02 

30,61 

3.4974E-03 

15,15 

3.4578E-04 

6,128 

1.1426E-03 

120,120 

1.2803E-03 

6,4,128 

2.6268E-04 

64,129 

9.4197E-05 

30,30 

4.7708E-05 

8,512 

1.0585E-04 

500,500 

1.1013E-04 

8,4,512 

3.5822E-05 

256,513 

8.8818E-16 

512,512 

9.9912E-09 


Table 2: The Relative Error for Example 1 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

8.8286E-02 

10,10 

4.1877E-02 

3,2,8 

8.8286E-02 

4,9 

1.8627E-01 

4,4 

1.4297E-03 

4,32 

3.1475E-03 

30,30 

3.5614E-03 

4,4,32 

2.2521E-03 

30,61 

7.7175E-04 

15,15 

7.6300E-05 

6,128 

2.5212E-04 

120,120 

12.825 IE-04 

6,4,128 

5.7963E-05 

64,129 

2.0786E-05 

30,30 

1.0527E-05 

8,512 

2.3356E-05 

500,500 

2.4302E-05 

8,4,512 

7.9046E-06 

256,513 

1.9599E-16 

512,512 

2.2047E-09 


Example 2 

r~ 1 sir? 2008 * 

f 2 1 _ SlU X _ fly 

(2007*+l) ( sin 2008 x+cos 2008 x ) 

The exact solution is ^ « 0.7853981633974483. Using eqn. (30). we calculated the absolute error, where x a is 
the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and exact 
solution are shown in Table 3. and using eqn. (31). we calculated the relative error and shown in Table 4. From the table as 
considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 2. obtained with 
above method. 
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Figure 2: Time Evolution of Given Methods for Example 2 


Table 3: The Absolute Error for Example 2 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine function 

Gauss Leg. Quad. 

(J,NN) 

e a 

(m,NN) 

e a 

(k,M,NN) 


(N,NN) 

e a 

(n,NN) 

e a 

2,8 

0.0000E+00 

10,10 

2.4758E-14 

3,2,8 

0.0000E+00 

4,9 

1.6641E-01 

4,4 

2.3899E-01 

4,32 

0.0000E+00 

30,30 

8.2157E-15 

4,4,32 

1.1102E-16 

30,61 

1.4478E-01 

15,15 

2.2177E-02 

6,128 

0.0000E+00 

120,120 

2.2204E-16 

6,4,128 

0.0000E+00 

64,129 

6.8996E-02 

30,30 

1.1467E-02 

8,512 

3.3307E-16 

500,500 

2.2204E-16 

8,4,512 

2.2204E-16 

256,513 

1.0589E-02 

512,512 

2.9933E-03 


Table 4: The Relative Error for Example 2.! 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

0.0000E+00 

10,10 

3.1523E-14 

3,2,8 

0.0000E+00 

4,9 

12.1188E-01 

4,4 

3.0429E-01 

4,32 

0.0000E+00 

30,30 

1.0460E-14 

4,4,32 

1.4136E-16 

30,61 

1.8434E-01 

15,15 

2.8237E-02 

6,128 

0.0000E+00 

120,120 

2.8272E-16 

6,4,128 

0.0000E+00 

64,129 

8.7848E-02 

30,30 

1.4601E-02 

8,512 

4.2407E-16 

500,500 

2.8272E-16 

8,4,512 

2.8272E-16 

256,513 

1.3482E-02 

512,512 

3.8112E-03 


4.2. Examples for Double Integration 
Example 3 

Jo 1 i 1 T^ dxd y 

The exact solution is Ym=i ~z ~ 1-644934066848226. Using eqn. (30). we calculated the absolute error, where 

x a is the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and exact 
solution are shown in Table 5. and using eqn. (31). we calculated the relative error and shown in Table 6. From the table as 
considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 3. obtained with 
above method. 
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Figure 3: Time Evolution of Given Methods for Example 3 


Table 5: The Absolute Error for Example 3 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e a 

(m,NN) 

e a 

(k,M,NN) 

e a 

(N,NN) 


(n,NN) 

e a 

2,8 

5.7222E-02 

10,10 

4.6446E-02 

3,2,8 

5.7222E-02 

4,9 

1.2976E-02 

4,4 

3.8940E-02 

4,32 

1.5186E-02 

30,30 

1.6173E-02 

4,4,32 

1.1349E-02 

30,61 

9.3798E-07 

15,15 

3.3315E-03 

6,128 

3.8718E-03 

120,120 

4.1279E-03 

6,4,128 

2.8602E-03 

64,129 

5.2915E-10 

30,30 

8.6129E-04 

8,512 

9.7397E-04 

500,500 

9.9729E-04 

8,4,512 

7.1653E-04 

125,251 

3.6193E-14 

512,512 

3.1996E-06 


Table 6: The Relative Error for Example 3 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

3.4787E-02 

10,10 

2.8236E-02 

3,2,8 

3.4787E-02 

4,9 

7.8883E-03 

4,4 

2.3673E-02 

4,32 

9.2320E-03 

30,30 

9.8321E-03 

4,4,32 

6.8992E-03 

30,61 

5.7022E-07 

15,15 

2.0253E-03 

6,128 

2.3538E-03 

120,120 

2.5094E-03 

6,4,128 

1.7388E-03 

64,129 

3.2169E-10 

30,30 

5.2360E-04 

8,512 

5.9210E-04 

500,500 

6.0628E-04 

8,4,512 

4.3560E-04 

125,251 

2.2003E-14 

512,512 

1.9451E-06 


Example 4 

J 2 f 4 x dxdy 

J 0 J x 2 7x 4 +y 2 * 

The exact solution is 4(V2 — 1) « 1.656854249492381. Using eqn. (30). we calculated the absolute error, 
where x a is the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and 
exact solution are shown in Table 7. and using eqn. (31). we calculated the relative error and shown in Table 8. From the table 
as considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 4. obtained 
with above method. 


Table 7: The Absolute Error for Example 4 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss 

Leg. Quad. 

(J,NN) 

e„ 

(m,NN) 

e a 

(k,M,NN) 

e„ 

(N,NN) 

e„ 

(n,NN) 

e a 

2,8 

9.6835E-03 

10,10 

5.9816E-03 

3,2,8 

9.6835E-03 

3,7 

8.4328E-03 

4,4 

7.8592E-04 

3,16 

2.1421E-03 

16,16 

2.1421E-03 

3,3,12 

6.4107E-04 

4,9 

9.7525E-04 

8,8 

1.4282E-04 

4,32 

4.5917E-04 

32,32 

4.5917E-04 

3,4,16 

4.7310E-04 

6,13 

4.3953E-04 

12,12 

2.9766E-05 

5,64 

9.5168E-05 

64,64 

9.5168E-05 

4,4,32 

1.4638E-04 

7,15 

3.3811E-04 

16,16 

9.4974E-06 
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Figure 4: Time Evolution of given Methods for Example 4 


Table 8: The Relative Error for Example 4 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

5.8445E-03 

10,10 

3.6102E-03 

3,2,8 

5.8445E-03 

3,7 

5.0896E-03 

4,4 

4.7435E-04 

3,16 

1.2929E-03 

16,16 

1.2929E-03 

3,3,12 

3.8692E-04 

4,9 

5.8861E-04 

8,8 

8.6197E-05 

4,32 

2.7713E-04 

32,32 

2.7713E-04 

3,4,16 

2.8554E-04 

6,13 

2.6528E-04 

12,12 

1.7965E-05 

5,64 

5.7439E-05 

64,64 

5.7439E-05 

4,4,32 

8.8347E-05 

7,15 

2.0407E-04 

16,16 

5.7322E-06 


4.3. Examples for Triple Integration 
Example 5 

Jo 1 Jo 1 Jo 1 TZ^- z dxd y dz 

The exact solution is Ym=i ~3 ~ 1-202056903159594. Using eqn. (30). we calculated the absolute error, where 

x a is the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and exact 
solution are shown in Table 9. and using eqn. (31). we calculated the relative error and shown in Table 10. From the table as 
considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 5. obtained 
with above method. 


Table 9: The Absolute Error for Example 5 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e a 

(m,NN) 

e a 

(k,M,NN) 

e a 

(N,NN) 

e a 

(n,NN) 

e a 

2,8 

4.0677E-03 

10,10 

2.8589E-03 

3,2,8 

4.0677E-03 

3,7 

1.0096E-02 

8,8 

8.8419E-05 

3,16 

1.3319E-03 

16,16 

1.3319E-03 

3,4,16 

4.5409E-04 

7,15 

1.0068E-03 

16,16 

6.3044E-06 

4,32 

4.1445E-04 

32,32 

4.1445E-04 

4,4,32 

1.1999E-04 

16,33 

1.821 IE-05 

32,32 

4.2024E-07 

5,64 

1.2435E-04 

64,64 

1.2435E-04 

4,8,64 

4.0084E-06 

32,65 

1.1358E-07 

64,64 

2.7112E-08 


Table 10: The Relative Error for Example 5 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

3.3839E-03 

10,10 

2.3783E-03 

3,2,8 

3.3839E-03 

3,7 

8.3988E-03 

8,8 

7.3556E-05 

3,16 

1.1081E-03 

16,16 

1.1081E-03 

3,4,16 

3.7776E-04 

7,15 

8.3758E-04 

16,16 

5.2447E-06 

4,32 

3.4479E-04 

32,32 

3.4479E-04 

4,4,32 

9.9824E-05 

16,33 

1.5150E-05 

32,32 

3.4960E-07 

5,64 

1.0345E-04 

64,64 

1.0345E-04 

4,8,64 

3.3346E-06 

32,65 

9.4487E-08 

64,64 

2.2555E-08 
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Example 6 

fo f/ 4 ~ y2 ~ z2 (x 2 +y 2 + Z 2 )dxdydz 

The exact solution is ^ « 10.053096491487338. Using eqn. (30). we calculated the absolute error, where x a is 

the approximate value of the integral with the above method and compared with Gauss Legendre quadrature and exact 
solution are shown in Table 11. and using eqn. (31). we calculated the relative error and shown in Table 12. From the table as 
considering NN be the number of nodes. An illustrate the computation performance of integration in Figure 6. obtained with 
above method. 



Table 11: The Absolute Error for Example 6 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 


(m,NN) 

e a 

(k,M,NN) 


(N,NN) 

e„ 

(n,NN) 


2,8 

9.7155E-02 

10,10 

6.8570E-02 

3,2,8 

9.7155E-02 

2,5 

5.7273E-02 

8,8 

5.1733E-03 

3,16 

3.2971E-02 

16,16 

3.2971E-02 

3,4,16 

1.5726E-02 

3,7 

4.1350E-02 

16,16 

6.9555E-04 

4,32 

1.1262E-02 

32,32 

1.1262E-02 

4,4,32 

5.4036E-03 

7,15 

6.0570E-03 

32,32 

9.0621E-05 

5,64 

3.8747E-03 

64,64 

3.8747E-03 

4,8,64 

7.7776E-04 

32,65 

6.5828E-07 

64,64 

1.1580E-05 
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Table 12: The Relative Error for Example 6 


Haar Wavelet 

BPF 

Chebyshev Wavelet 

Sine Function 

Gauss Leg. Quad. 

(J,NN) 

e r 

(m,NN) 

e r 

(k,M,NN) 

e r 

(N,NN) 

e r 

(n,NN) 

e r 

2,8 

9.6642E-03 

10,10 

6.8208E-03 

3,2,8 

9.6642E-03 

2,5 

5.6971E-03 

8,8 

5.1459E-04 

3,16 

3.2797E-03 

16,16 

3.2797E-03 

3,4,16 

1.5643E-03 

3,7 

4.1131E-03 

16,16 

6.9188E-05 

4,32 

1.1203E-03 

32,32 

1.1203E-03 

4,4,32 

5.3750E-04 

7,15 

6.0250E-04 

32,32 

9.0143E-06 

5,64 

3.8543E-04 

64,64 

3.8543E-04 

4,8,64 

7.7366E-05 

32,65 

6.5481E-08 

64,64 

1.1519E-06 


5. RESULTS AND DISCUSSIONS 

Maximum absolute and relative errors obtained by the numerical multiple integrations based on the given methods 
for certain problems are shown in Tables 1 to 12. Comparison of maximum absolute and relative errors with Gaussian 
Legendre quadrature method and exact solution determines the accuracy of given methods. Estimates the computation 
performance for the certain problems is shown in figure 1 to 6. Comparison of time evolutions with Gaussian Legendre 
quadrature method. It is clearly seen that numerical multiple integrations using above methods have produced improved 
results ,as well as computation performance, are compared with Gaussian Legendre quadrature method for numerical 
multiple integrations. 

6. CONCLUSIONS 

In this paper, Haar wavelet, Block-pulse function, Chebyshev wavelet and Sine functions were applied for 
numerical multiple integrations (double, triple integrals) with variable and constant limits. Compared between these methods 
with Gaussian Legendre quadrature method. The exact solution clearly shows that the above methods and given a much 
better (anyone of Haar wavelet, Block-pulse function, Chebyshev wavelet and Sine functions and it depends on the problem) 
results as well as computation performance than Gaussian Legendre quadrature method. 
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